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ABSTRACT 


In  this  report,  an  evaluation  of  the  Watson- Watt  and  Butler  matrix  approaches 
for  tactical  wideband  radio  direction  finding  applications  is  described.  This  evaluation 
was  carried  out  using  theoretical  derivations  and  computer  simulations,  concentrating  on 
the  effects  of  various  error  mechanisms  including  internal  noise,  terrestrial  and  isotropic 
noise,  cochannel  interference,  and  multipath.  A  maximum  likelihood  (ML)  approach  was 
also  derived  and  evaluated  for  comparative  purposes  to  determine  the  accuracy  trade-off 
between  using  an  A’-channel  approach  versus  the  three  channel  Watson- Watt  approach 
and  the  two  channel  Butler  matrix  approach.  The  results  indicate  that  the  Watson- 
Watt  approach  is  superior  to  the  Butler  matrix  approach  for  current  applications.  In  the 
future,  this  choice  may  need  to  be  reassessed  as  lower  cost  receivers  and  faster  processing 
equipment  make  more  sophisticated  approaches  such  as  ML  attractive  alternatives. 

RESUME 

Le  present  resume  decrit  revaluation  des  approches  fondees  sur  Watson- Watt  et 
la  matrice  de  Butler  pour  les  applications  radiogoniometriques  tactiques  a  large  bande. 
L’evaluation  a  fait  appel  a  des  calculs  theoriques  et  a  des  simulations  par  ordinateur  por- 
tant  sur  les  effets  des  divers  m^canismes  d’erreur,  dont  le  bruit  interne,  le  bruit  terrestre  et 
isotopique,  I’interference  entre  voies  et  la  propagation  par  trajets  multiples.  L’evaluation 
a  egalement  permis  de  deriver  et  d’evaluer  une  approche  basee  sur  la  probabilite  maximale 
(PM)  pour  comparer  les  precisions  angulaires  obtenues  en  utilisant  N  canaux,  les  trois 
canaux  de  Watson- Watt  et  une  matrice  Butler  a  deux  canaux.  Les  resultats  montrent 
que  I’approche  de  Watson- Watt  est  superieure  a  la  matrice  de  Butler  pour  les  applica¬ 
tions  courantes.  Dans  I’avenir,  on  devra  peut-etre  reevaluer  ce  choix,  car  les  recepteurs 
peu  couteux  et  le  materiel  de  traitement  accelere  rendent  plus  interessantes  les  approches 
poussees  comme  celles  de  la  PM. 


EXECUTIVE  SUMMARY 


Advances  in  technology  have  seen  the  development  of  improved  communica¬ 
tions  receivers  with  larger  and  larger  bandwidths.  Coupled  with  advances  in  digital  and 
computing  technology,  these  new  wideband  receivers  allow  surveillance  of  spectral  bands 
containing  many  communications  signals  with  100%  detection  of  all  transmissions  in  that 
band.  The  development  of  tactical  wideband  surveillance  systems  has  also  led  to  an  in¬ 
terest  in  wideband  radio  direction  finding  (DF)  with  the  ultimate  goal  of  locating  the 
source  of  transmission.  To  keep  equipment  size  and  cost  to  a  minimum,  DF  approaches 
have  been  proposed  which  minimize  the  number  of  receiver  channels  by  performing  analog 
processing  at  the  front  end. 

Two  promising  implementations  for  wideband  DF  are  the  Watson- Watt  and  the 
Butler  matrix  approaches.  Although  philosophically  different,  they  have  similar  hardware 
and  processing  requirements.  To  evaluate  and  compare  these  two  approaches  theoretical 
derivations  and  computer  simulations  were  carried  out  where  the  effects  of  noise,  cochannel 
interference  and  multipath  were  present.  Additionally,  the  performances  of  these  two 
approaches  were  also  compared  with  a  more  optimum  maximum  likelihood  (ML)  approach 
to  determine  the  accuracy  penalty  incurred  by  reducing  the  number  of  receiver  channels 
(the  ML  approach  requires  one  receiver  for  every  antenna). 

For  the  theoretical  evaluations  of  the  Watson- Watt  and  Butler  matrix  approaches, 
four  and  eight-element  circular  antenna  arrays  were  considered.  It  was  found  that  when 
adjusting  the  size  of  these  arrays,  two  critical  antenna  array  radii  exist  where  the  gains 
of  one  or  more  of  the  receiver  channels  becomes  zero  leading  to  large  bearing  errors 
when  noise  or  other  types  of  interference  are  present.  The  smallest  radius  (approximately 
one-third  the  signal  wavelength)  at  which  this  happens  corresponds  to  zero  gain  in  the 
reference  channel  (a  channel  common  to  both  approaches  and  used  for  phase  reference 
purposes).  The  problem  can  be  overcome  by  using  a  central  antenna  for  the  reference 
instead  of  deriving  it  from  the  antennas  in  the  circular  array,  although  this  could  be  dif¬ 
ficult  to  do  if  a  supporting  mast  occupies  the  central  position.  The  second  critical  radius 
(approximately  one-half  the  signal  wavelength  or  slightly  more)  corresponds  to  the  gain 
in  the  other  channels  going  to  zero,  and  there  is  no  simple  solution.  Hence  the  maximum 
size  is  limited  by  this  second  critical  radius.  For  a  fixed  array  size,  the  upper  frequency 
is  limited  accordingly.  The  ML  approach  is  also  limited  by  the  same  critical  radii  for 
the  four-element  array  but  not  for  the  eight-element  array  which  is  capable  of  working 
up  to  an  array  radius  of  slightly  less  than  three-quarters  of  the  signal  wavelength  for  an 
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eight-element  antenna  array  with  a  correspondingly  higher  upper  frequency  limit. 

In  practical  applications,  external  effects  such  as  terrestrial  and  cosmic  noise, 
cochannel  interference,  and  multipath  will  be  the  dominant  source  of  error.  Under  these 
conditions,  the  accuracies  of  both  the  Watson- Watt  and  Butler  matrix  approaches  were 
determined  to  be  relatively  constant  as  the  array  size  was  increased  (or  equivalently  as 
the  frequency  was  increased)  except  close  to  the  critical  values  (discussed  in  the  previous 
paragraph)  where  accuracy  rapidly  degraded.  There  was  relatively  little  improvement  in 
either  of  the  two  approaches  when  an  eight-element  antenna  array  was  used  compared  to 
a  four-element  array.  Quantitatively,  the  Watson- Watt  approach  was  found  to  be  twice  as 
accurate  (i.e.  one  half  the  RMS  bearing  errors)  as  the  Butler  matrix  approach.  In  com¬ 
parison,  the  ML  approach  was  significantly  better  with  the  best  accuracy  achieved  using 
the  eight-element  array  (upwards  of  five  times  better  than  the  Watson- Watt  approach). 

Generally  the  results  show  that  neither  the  Watson- Watt  nor  Butler  matrix  ap¬ 
proaches  achieve  optimum  accuracy  under  practical  conditions  and  that  utilizing  more 
channels  provides  better  accuracy.  However,  for  a  system  where  other  constraints  be¬ 
sides  accuracy  are  important,  such  as  low  cost  and  realtime  capability,  and  given  the 
current  economic  climate  and  technological  advances,  the  Watson- Watt  approach  is  the 
best  choice.  The  ML  approach,  although  generally  having  better  accuracy  and  a  greater 
frequency  range  of  operation,  is  considerably  more  expensive  due  to  the  greater  number  of 
receivers  and  the  processing  requirements.  In  the  future,  however,  this  choice  may  need 
to  be  reassessed  as  lower  cost  receivers  and  faster  processing  equipment  make  the  ML 
approach  a  more  attractive  alternative. 
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1.0  INTRODUCTION 


Advances  in  technology  have  seen  the  development  of  improved  communications  re¬ 
ceivers  with  larger  and  larger  bandwidths.  Coupled  with  advances  in  digital  and  comput¬ 
ing  technology,  these  new  wideband  receivers  allow  surveillance  of  spectral  bands  contain¬ 
ing  many  communications  signals  with  100%  detection  of  all  transmissions  in  that  band 
-  something  that  was  not  possible  previously.  Once  detected,  these  signals  may  then  be 
demodulated  and  analyzed  as  necessary. 

The  development  of  tactical  wideband  surveillance  systems  has  also  led  to>an  interest 
in  wideband  radio  direction  finding  (DF)  with  the  ultimate  goal  of  locating  the  source 
of  transmission.  To  keep  equipment  size  and  cost  to  a  minimum,  DF  approaches  have 
been  proposed  which  minimize  the  number  of  receiver  channels  by  performing  analog 
processing  at  the  front  end.  The  front  end  processing,  however,  leads  to  compromises  in 
performance  since  there  is  less  flexibility  in  terms  of  how  the  received  antenna  signals  are 
processed. 

Two  approaches  to  wideband  DF,  which  appear  to  have  the  desired  qualities,  are  the 
Watson- Watt  approach  and  the  Butler  matrix  approach.  Although  philosophically  dif¬ 
ferent  approaches,  their  hardware  implementation  is  similar  and  one  might  expect  that 
the  performance  of  both  approaches  would  also  be  similar.  However,  subtle  differences 
in  approaches  can  often  lead  to  unexpected  and  undesired  consequences  in  terms  of  per¬ 
formance  which  is  not  something  one  wishes  to  discover  after  the  equipment  has  been 
built.  To  address  these  concerns,  and  identify  the  more  promising  approach,  an  evalua¬ 
tion  of  these  two  approaches  was  carried  out  through  theoretical  derivations  and  computer 
simulations.  This  evaluation  included  the  effects  of  noise,  cochannel  interference  and  mul¬ 
tipath.  Additionally,  the  performance  of  the  Watson- Watt  and  Butler  matrix  approaches 
was  compared  with  a  more  optimum  maximum  likelihood  approach  to  highlight  the  per¬ 
formance  trade-off  being  made. 

This  introduction  is  followed  by  a  theoretical  derivation  of  the  Watson- Watt  approach, 
the  Butler  matrix  approach,  and  the  maximum  likelihood  (ML)  approach  in  Section  2. 
This  is,  in  turn,  followed  by  an  evaluation  of  the  theoretical  and  simulated  susceptibility  of 
the  proposed  wideband  approaches  to  internal  noise  and  external  interference  (including 
noise,  cochannel  interference,  and  multipath)  in  Section  3.  The  simulated  susceptibil¬ 
ity  of  the  ML  approach  is  also  included  for  comparison.  Finally,  the  conclusions  and 
recommendations  are  presented  in  Section  5. 


2.0  DF  THEORY 


For  the  purposes  of  this  report  only  four  and  eight-element  equally  spaced  circular 
antenna  arrays  are  considered.  This  does  not  mean  that  the  Watson- Watt  or  Butler 
matrix  approaches  are  restricted  to  only  these  two  antenna  arrays,  rather  it  provides 
a  convenient  basis  for  comparison  and  these  are  commonly  used  arrays  for  tactical  DF 
applications. 

Since  none  of  the  DF  approaches  discussed  are  capable  of  dealing  with  more  than 
one  signal  at  a  time,  it  is  assumed  that  the  receiver  outputs  are  digitized  then  processed 
using  narrowband  digital  filters  or  FFT  processing  to  separated  the  various  signals  before 
DF  processing.  Regardless  of  how  a  target  signal  has  been  isolated,  the  input  to  the 
DF  processor  is  assumed  to  be  a  single  complex  sample  from  each  channel  at  a  sampling 
rate  considerably  lower  than  the  sampling  rate  of  the  receivers  (which  is  related  to  the 
difference  in  the  receiver  bandwidth  and  the  signal  bandwidth). 

Averaging  of  the  estimated  signal  bearing  may  also  be  done  to  produce  a  more  stable 
estimate,  however,  this  is  not  discussed  in  this  report.  Rather  it  is  assumed  that  if  the 
bearing  errors  have  a  random  temporal  component  the  error  distribution  will  be  Gaussian 
and  averaging  will  then  reduce  the  temporal  RMS  error  by  a  factor  of  1/\/K  where  K  is 
the  number  of  bearings  averaged.  Note  that  in  some  instances,  such  as  with  multipath  at 
the  receiver  site,  the  resultant  error  does  not  change  over  time  so  that  averaging  provides 
no  improvement. 

To  maintain  consistency  in  the  derivations,  the  positions  (x^,  y„)  for  n  =  0, 1, ...,  A"  —  1 
of  the  N  antennas  used  in  the  circular  arrays  are  defined  according  to 

x„  =  rcos^^H-i?^  (1) 

yn  =  ^sin^^-l-1?^  (2) 

where  the  array  has  radius  r,  the  origin  is  the  center  of  the  array,  and  -d  gives  the  angle 
of  the  radial  running  through  antenna  0  measured  counter-clockwise  from  the  X-axis  (the 
array  orientation).  Figures  1  and  2  show  the  antenna  positions  graphically.  The  response 
of  antenna  n  for  a  signal  arriving  from  the  direction  <f)  (measured  in  the  clockwise  direction 
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from  the  y-axis)  is  then  given  by 


Xfi 


se^^0^n^n(i>+y^cos^) 


=  sei¥-i"(*+*r+<')  (3) 

where  s  is  the  time  dependent  complex  amplitude  of  the  signal  measured  at  the  center  of 
the  antenna  array,  and  A  is  the  signal  wavelength. 


Figure  1:  Geometry  of  a  four-element  circular  array  showing  antenna  positions  (large 
dots)  and  orientation  =  15°,  as  well  as  signal  direction  <l>  =  45°  (where  the  signal  is 
represented  by  the  large  arrow). 


Before  continuing  on  to  the  detailed  derivation  of  either  the  Watson- Watt  approach  or 
the  Butler  matrix  approach,  it  is  useful  to  consider  some  of  the  common  features.  Both 
approaches  involve  estimation  of  the  signal  direction  using  the  equation 

^  =  arctan 

where  A  and  B  are  real-valued  and  the  four  quadrant  version  of  the  arctangent  function 
is  used  (i.e.,  the  signs  of  A  and  B  are  utilized  to  resolve  the  180°  ambiguity  between  the 
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X-coordinate  (wavelengths) 


Figure  2:  Geometry  of  a  eight-element  circular  array  showing  antenna  positions  (large 
dots)  and  orientation  d  =  22.5°,  as  well  as  signal  direction  =  45°  (where  the  signal  is 
represented  by  the  large  arrow). 


quantities  4  and  The  four  quadrant  arctangent  function  also  has  the  property 


^  =  arctan 


(i) 


where  c  is  a  positive  real  value. 


In  determining  the  signs  of  the  quantities  A  and  B,  a  reference  channel  is  required. 
This  channel  is  usually  generated  by  summing  the  outputs  from  all  the  antennas  in  the 
array.  The  result  is  given  by 

N-l 

“^ref  ~  ^  ^ 

n=0 

N/2-1 

=  ^  ^  Xfi  +  Xji^pf/2 

n=0 

N/2-1 


=  sin(0+2^+7r+i?) 

n=0 

/27rr  .  27rn  ..N 

=  2s  ^  cos  sm((^ -I- — -fi?)  j 
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=  2s/3 


(6) 


where  is  a  real  scaling  factor.  The  value  of  0  changes  significantly  with  r,  the  radius  of 
the  antenna  array,  as  shown  in  Figures  3  and  4,  and  is  relatively  independent  of  the  signal 
direction  ^  except  for  values  of  r  >  0.3A  in  the  four  element  case.  Prom  these  figures,  the 
minimum  value  of  /3  occurs  for  ^  +  ??  =  (2A:  +  l)7r/N  for  k  =  0, ...,  iV  -  1  and  is  given  by 


Pmin  —  S 


'  f  TrrV2\ 

2cos(^— j 


-  /27rr  37r\  . 

2  cos  cos  —  1  +  2  cos 


/27rr  TT' 


for  N  =  4 


for  N  =  8 


(7) 


Similarly,  the  maximum  value  of  ^  occurs  for  ^  =  2kTr/N  for  k  =  0, ...,  N  -1  and  is 

given  by 

''27rr\ 


Pmax  —  A 


1  +  cos 


A  J 


for  N  =  4 


1  +  2  cos  j  +  cos  for  N  =  8 


(8) 


Alternatively,  for  the  Watson- Watt  approach,  the  reference  can  be  obtained  from  the 
output  of  a  single  antenna  placed  at  the  center  of  the  array.  The  reference  output  in  this 
case  is  simply 

Vref  —  S  (9) 

The  advantage  of  using  a  single  antenna  is  that  it  does  not  suffer  from  the  zero  problem 
associated  with  the  summed  output.  The  disadvantages  are  that  an  extra  antenna  and 
associated  hardware  are  required  and  the  central  position  of  the  array  will  typically  be 
occupied  by  the  mast  used  to  support  the  antenna  array.  Hence,  when  discussing  the 
reference  voltage,  the  form  given  by  (6)  is  assumed  unless  otherwise  stated. 

In  Figures  3  and  4,  goes  to  zero  then  becomes  negative  for  larger  values  of  r. 
Since  this  is  an  undesirable  condition  (particularly  ^  =  0),  a  suitable  limit  on  r  can  be 
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Figure  3:  Value  of  as  a  function  of  both  the  signal  bearing  and  the  radius  of  a  four- 
element  circular  antenna  array  (i9  =  0"). 
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Figure  4:  Value  of  jS  as  a  function  of  both  the  signal  bearing  and  the  radius  of  an 
eight-element  circular  antenna  array  (i?  =  0°). 
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determined  by  requiring  that  Prnm  >  0  which  leads  to 


'^maxo 


A 

2^/2 

A 

,  2-\/2cos| 


for  N  =  A 


for  iV  =  8 


(10) 


These  limitations  on  r  only  apply  when  the  reference  voltage  is  generated  using  the 
summed  antenna  outputs.  For  the  central  reference  antenna  approach  there  is  no  such 
limitation. 


It  is  also  useful  to  evaluate  the  term  Xn—Xn+N/2  which  figures  prominently  in  generating 
the  channel  outputs  for  both  approaches.  Hence 

Xn-Xn+N/2  =  56^'^ 

=  -  se-^^ 


.  /27rr  .  , ,  27rn  ^  , 

=  j2ssin(^— sm(^+  — +  79)j  (11) 

Finally,  some  comments  on  mutual  coupling  effects  are  in  order.  Given  the  small 
spacings  between  antennas  used  in  the  antenna  arrays  discussed  in  this  study  (less  than 
a  wavelength),  mutual  coupling  effects  would  be  expected  to  figure  prominently  in  the 
results.  However  it  was  found  through  simulation  that,  although  mutual  coupling  does 
significantly  affect  channel  gains  and  phases,  the  effects  could  be  minimized  through 
channel  calibration  with  the  result  that  performance  was  not  appreciably  different  than 
when  mutual  coupling  was  absent.  For  example,  assuming  half-wave  dipole  elements 
and  using  the  approach  discussed  in  [1]  to  model  mutual  coupling,  the  summed  reference 
voltage  amplitude  was  reduced  by  a  factor  ranging  from  N  (for  r  ^  0)  to  1  (for  r  Tmaxo)- 
The  value  of  Tmaxo  was  not  significantly  changed,  however,  and  the  phase  offset  was  easily 
corrected,  minimizing  the  degradation  of  the  reference  voltage  (the  reduction  in  amplitude 
was  minimal  for  values  of  r  where  bearing  estimation  could  have  been  compromised). 
Hence  to  avoid  unnecessarily  complicating  the  derivation  and  evaluation  of  the  direction 
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finding  approaches  discussed  in  this  report,  mutual  coupling  effects  were  ignored. 


2.1  The  Watson- Watt  Approach 

The  Watson- Watt  approach  uses  two  orthogonal  directional  antennas  with  the  gain 
pattern  of  one  antenna  given  by 

=  sin(^  (12) 

and  the  gain  of  the  second  antenna  given  by 


9ns{(I>)  =  cos  (j) 


(13) 


where  in  this  discussion  the  gains  are  real  valued  but  may  be  positive  or  negative.  If  these 
two  gain  values  are  measured  for  a  given  signal,  as  shown  in  Figure  5,  then  the  signal 
bearing  can  be  estimated  using 


(j)  =  arctan 


(14) 


Figure  5:  Antenna  patterns  for  an  eight-element  circular  array  implementing  the  Watson- 
Watt  approach. 


The  next  two  sections  discuss  the  implementation  of  the  Watson- Watt  approach  using 
four  and  eight-element  circular  arrays. 


2.1.1  Four  Element  Array 

Directional  antennas  with  the  sine  and  cosine  patterns  for  an  AT  =  4  element  circular 
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Figure  6:  Adcock  antenna  setup  for  a  four-element  array  Watson- Watt  system. 


array  can  be  approximated  by  taking  the  difference  in  outputs  from  opposite  pairs  of  an¬ 
tennas  as  shown  in  Figure  6  and  positioning  antenna  0  on  the  X-axis  (the  East- West  line). 
Using  (11)  and  the  fact  that  i9  =  0°,  the  sine  response  gain  pattern  can  be  approximated 
by  taking  the  difference  between  antennas  0  and  2  giving 


vew 


Xq  -  X2 


(15) 


where  vew  is  the  resultant  complex  voltage.  The  cosine  response  gain  pattern  can  be 
generated  in  a  similar  fashion  by  taking  the  difference  between  antennas  1  and  3  giving 


vns 


Xi  -  Xz 


(16) 


where  v^s  is  the  resultant  complex  voltage. 


Inspection  of  the  above  expressions  reveals  that  for  particular  combinations  of  the 
bearing  0  and  array  radius  r,  both  voltages  will  be  zero  (e.g.  =  (0°,A/2)  or 

(45®,  A/ v^)).  This  condition  can  be  avoided  if  the  array  radius  is  kept  small  enough 
that  r  <  rmaxi  where 


^tnaxi 


A 

2 


(17) 


The  condition  r  <  r^axi  then  ensures  that  the  directional  signal  voltages  will  be  nonzero 
regardless  of  the  signal  bearing,  although  for  some  signal  bearings  it  will  be  possible  to 
exceed  this  value.  However,  unless  the  reference  voltage  is  derived  from  the  output  of  a 
central  antenna  instead  of  using  summed  antenna  outputs  as  described  in  Section  2.0,  the 
actual  upper  limit  on  the  array  size  will  be  r  <  r^axo  (where  Tmaxo  is  given  in  (10))  since 

'^maxo  ^  ‘^maxi’ 
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Although  (15)  and  (16)  are  not  exactly  the  desired  response,  if  the  antenna  array  size 
is  small  enough,  then  the  above  expressions  reduce  to  the  more  desirable  forms 


and 


respectively. 


QiTtt 

vew  ~  3‘^s——  sin  (j)  for  r  A/27r 
A 


27r7’ 

vns  ^  j^s——  cos  <}>  for  r  ■<  X/2'k 
A 


(18) 


(19) 


To  convert  the  complex  voltages  to  more  usable  forms  appropriate  for  a  four  quadrant 
arctangent  function,  since  in  most  practical  applications  only  relative  phases  can  be  mea¬ 
sured,  then  the  two  voltages  are  phase  shifted  -90®  and  divided  by  the  reference  voltage 
Vrej  yielding 


and 


_  .Vew 

1  . 

=  -sm 

r^sin<?i) 

(20) 

'^ref 

P 

\  A  / 

-  _ 

■-  4  sin  ( 

'27rr 

— —  COS  o 

(21) 

'^ref 

P  V 

.  A  J 

Provided  that  r  <  A/2-\/2  then  the  term  P  will  always  be  positive  and  real,  and  g'Ew{4>) 
and  g'ifs{4>)  can  be  used  to  replace  gEw{4>)  and  gNs{<l>)  respectively,  in  (14).  Hence 


0  =  arctan 


^sin  ^^sin^^ 
^sin  cos  (/>) 


(22) 


The  approximations  (15)  and  (16)  used  for  the  cosine  patterns  lead  to  bias  errors  which 
are  a  function  of  both  signal  direction  and  antenna  array  size  as  shown  in  Figure  7.  As 
noted  already,  as  the  antenna  array  size  shrinks,  the  approximations  improve  and  the  bias 
error  is  reduced.  Since  the  bias  errors  are  constant  and  repeatable,  a  look-up  table  can 
be  used  to  correct  the  results  with  relatively  little  processing  overhead. 
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Figure  7:  Bias  errors  of  a  Watson- Watt  system  using  a  four-element  array  showing  errors 
as  a  function  of  the  radius  of  the  antenna  array  and  the  signal  bearing. 

2.1.2  Eight  Element  Array 


Figure  8:  Adcock  antenna  setup  for  an  eight-element  array  Watson- Watt  system. 


Directional  antennas  with  the  sine  and  cosine  patterns  for  an  AT  =  8  element  circular 
array  can  be  approximated  by  adding  adjacent  pairs  of  antennas  and  taking  the  difference 
in  outputs  between  these  antennas  and  their  opposite  counterparts.  To  properly  align 
these  patterns  with  the  X  and  Y  axes,  the  array  is  oriented  so  =  — tt/S  (this  aligns  the 
long  axis  of  rectangle  formed  using  antennas  0, 1,  4,  and  5  with  the  X-axis).  The  resultant 
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voltage  outputs  in  this  case  are  given  by 


VeW  —  (^0  +  ^l)  ~  (^4  + 


=  32s  sin  sin(</>  -  +  j2s  sin  sin(.^  + 


cos(^)  sin  (t>  for  r  *C  X/2-k 

A  8 


VnS  —  (^2  +  ^3)  ~  (^6  +  ^7) 

=  j2s  sin  cos{<f>  -  +  ;2s  sin  cos(0  +  (25) 

~  js-^cos(^)cos<l)  for  r  «C  A/27r  (26) 

A  8 

where  =  tt/S  was  used  to  adjust  the  array  so  that  the  long  side  of  the  rectangle  formed 
using  antennas  0,  1,  4,  and  5  was  aligned  with  the  X-axis  (which  also  aligns  the  resultant 
sine  and  cosine  voltage  patterns  in  the  proper  directions). 

As  in  the  four-element  case,  certain  combinations  of  the  signal  bearing  (f)  and  the  an¬ 
tenna  array  radius  result  both  channel  voltages  becoming  zero  (e.g.  (0,  r)  =  (0°,  A/2  cos  |) 
or  (45°,  A/V2)).  This  condition  can  be  avoided  if  the  array  radius  is  kept  small  enough 
that  r  <  Tmaxi  where 

'  ""  2^ 

This  upper  limit  on  the  array  radius  is  slightly  larger  than  the  upper  limit  given  for  the 
four-element  array  in  (17)  and  applies  for  the  case  when  a  reference  antenna  is  used  at 
the  center  of  the  array.  For  the  reference  voltage  generated  using  the  summed  antenna 
output  the  upper  limit  given  by  (10)  is  smaller  and  therefore  applies. 

Shifting  the  directional  antenna  voltages  -90°  and  dividing  by  the  reference  voltage 
to  get  the  relevant  gain  values  then 
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and 


9'Nsi<f>)  =  =  ^sin^^cos(<A-|)) +isin^^cos((^  +  ^))  (29) 


Finally,  the  bearing  estimate  is  determined  using 

^=arctan('^2#') 

\9ns{<p)  ) 

In  terms  of  (28)  and  (29)  the  above  expression  is  given  by 

^  nrctan  (  "  f ))  +  ^“  (t  +  f )) 

ysin  cos(0  -  |))  +  sin  cos(^  +  f )) 


(30) 


(31) 


Figure  9:  Bias  error  of  a  Watson- Watt  system  using  an  eight-element  array  as  a  function 
of  the  radius  of  the  antenna  array  and  signal  bearing. 

As  in  the  four  element  case,  the  approximations  (23)  and  (25)  used  for  the  cosine 
patterns  lead  to  significant  bias  errors  in  the  bearing  estimates  except  when  the  array 
size  is  very  small  as  shown  in  Figure  9.  However  the  magnitude  of  the  errors  in  this  case 
are  about  one-third  those  for  the  four-element  system.  Since  the  bias  errors  are  constant 
and  repeatable,  a  look-up  table  can  be  used  to  correct  the  results  with  relatively  little 
processing  overhead. 
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2.2  The  Butler  Matrix  Approach 


The  Butler  matrix  approach  takes  advantage  of  the  relationship  between  the  spatial 
response  measured  around  a  circular  array  and  the  signal  direction.  For  example,  the  spa¬ 
tial  phases  measured  around  a  one  wavelength  diameter  circle  for  various  signal  bearings 
are  shown  in  Figure  10  (the  corresponding  amplitude  values  are  not  shown  since  they  are 
constant  regardless  of  signal  direction  or  angular  position).  As  can  be  seen,  the  shape  of 
the  spatial  response  is  unaffected  by  the  signal  direction,  so  that  the  signal  direction  can 
be  related  directly  to  the  phase  of  this  response  waveform. 


Figure  10:  Spatial  phase  response  measured  around  a  one  wavelength  diameter  circle 
for  signal  directions  (j)  =  0°,  60°,  120°,  and  180°.  Phase  measurements  were  made  relative 
to  the  center  of  the  circle. 

Using  a  circular  array  to  make  discrete  measurements,  one  method  of  determining  the 
spatial  phase  of  the  response  waveform  is  to  use  a  Butler  matrix  which  performs  a  discrete 
spatial  Fourier  transform  on  the  data.  Treating  the  antenna  outputs  as  N  points  in  a 
linear  spatial  sequence,  the  transform  produces  the  sequence 

(32) 

where  /o,  /i, ...,  Sn-i  S'l’e  the  Fourier  frequency  coefficients.  These  outputs  are  also  called 
the  mode  voltages  where  the  mode  0  voltage  is  /o,  the  mode  1  voltage  is  /i,  and  so  on. 
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The  coefficients  are  calculated  using 


2nik 


N-1 

fk=Y^ 

i=0 


(33) 


Given  that  one  complete  waveform  requires  all  N  antennas,  and  ignoring  the  arbitrary 
phase  term  introduced  by  s,  the  signal  modulation,  then  the  Fourier  coefficient  fi  contains 
the  relevant  spatial  phase  information.  Since  the  modulation  term  s  affects  each  Fourier 
term  in  the  same  way,  its  effect  can  be  removed  by  dividing  through  by  /q.  The  spatial 
phase  will  then  be  the  phase  of  fi/fo.  Since  the  phase  will  change  exactly  as  the  signal 
direction  changes,  then  for  a  suitable  choice  of  the  array  orientation  (i9  =  0),  the  phase 
angle  and  the  estimated  signal  bearing  will  be  the  same.  Hence  the  estimated  signal 
bearing  will  be  given  by 


(f)  =  arctan 


Amag{/i//ii}\ 
I  real{/,//„}  } 


(34) 


2.2.1  Four  Element  Array 

From  (33),  the  mode  0  voltage  fo  is  the  sum  of  the  outputs,  or  Vref,  and  is  calculated 
according  to  (6)  where  =  0.  The  Fourier  coefficient  for  fi  is  given  by 

fi  =  XQ+  xie~^^  +  +  xze~^^ 


=  Xq-  jXi  -X2+  jXz 


=  {xq  -  X2)  -  -  Xz) 


Expanding  the  two  bracketed  terms  using  (11)  for  n  =  0  and  n  =  1  then 


h 


0) 


+  2s  sin 


+  2s  sin 


Dividing  this  result  by  /o  =  Vrej,  where  is  given  by  (6),  then 


/1//0 


+  j-sin 


(35) 


(36) 
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Taking  (37)  and  plugging  it  into  (34)  results  in 


sin  (to  dint)  \ 

sm(tocos^)j  ^  ' 

which  is  identical  to  (22),  the  expression  developed  for  the  Watson- Watt  approach  using 
four  antennas.  Hence  the  bias  error  due  to  approximating  a  continuous  circle  by  four 
discrete  points  (i.e.  using  an  array  of  four  elements)  is  given  by  Figure  7. 


^  =  arctan 


Much  like  the  problem  of  channel  voltages  going  to  zero  in  the  Watson- Watt  approach, 
the  mode  voltage  fo  and  fi  also  go  to  zero  for  particular  combinations  of  and  r.  The 
zero  condition  for  fo  =  Vref  was  previously  discussed  in  Section  2.0  while  the  zero  condition 
for  fi  is  identical  to  the  four-element  Watson- Watt  approach  discussed  in  Section  2.0  and 
the  same  limitations  on  array  size  also  apply.  Additionally,  although  a  central  reference 
antenna  is  not  normally  associated  with  the  Butler  matrix  approach,  this  does  not  preclude 
it’s  use  to  increase  the  usable  size  of  the  antenna  array  to  r  <  rmax2  where  rmax2  =  f'maxi  • 


2.2.2  Eight  Element  Array 

The  Fourier  coefficient  for  fi  for  an  eight  element  array  is  given  by 

fi  =  xo  +  xie~^^  -I-  X2e~^2  -I-  Xze~^^  -I-  Xie~^‘"  -}-  +  Xoe~^^  -1-  X7e~^^ 


,  1  —  i  .  1+j  1  “  i  ,  .  ,  1  +  i 

Xq  +  -X4-  -^^5  +  JX6  +  -^^7 

(xo  -  X4)  +  “^(®1  ~  ^5)  -  J(X2  -  Xe)  -  ”^(^3  -  X7) 


(xo  -  X4)  +  ^(^1  -  ^5)  ^7) 

-  ^5)  -  J(x2  -  xe)  -  -  X7) 


Expanding  the  bracketed  terms  using  (11)  for  n  =  0, 1, 2, 3  and  then  dividing  through  by 


47rr  ...  A 

«  -TTr  sm  ffl  for  r  -C  — 

Xp  27r 


(44) 


Substituting  (41)  and  (43)  into  (34)  leads  to  the  result 

f ))  +  V2  sin  (y  sin  <f>)  +  sin  cos{(f)  +  f ))  \ 

f ))  +  \/2sin  cos  <f>^  +  sin  sin(^  +  f ))  ) 

which  is  obviously  quite  different  than  the  corresponding  result  (31)  for  the  Watson- Watt 
system.  The  bias  error  (due  to  discrete  spatial  sampling)  is  almost  insignificant  in  this 
case  as  shown  in  Figure  11.  Only  the  result  for  r  =  0.3 A  was  shown  since  for  smaller 
values  of  r  the  biasing  is  even  smaller. 
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Figure  11:  Bias  error  of  a  Butler  matrix  system  using  an  eight-element  array  as  a  function 
of  the  signal  bearing. 

The  problem  of  fo  and  fi  going  to  zero  for  certain  combinations  of  bearings  and  array 
sizes  also  affects  the  eight-element  Butler  matrix  approach.  The  limitation  on  the  array 
size  based  on  |/o|  =  \vref\  >  0  is  given  by  (10),  while  the  limitation  on  the  array  size  based 
on  |/i|  >0  is  given  by 

rmax2  =  0.60564751A  (46) 

Since  the  above  result  is  larger  than  given  by  (10),  the  above  limit  is  superseded  by 
(10).  Although  it  would  require  a  modification  to  the  basic  Butler  matrix  approach,  the 
above  limit  could  be  reached  by  replacing  the  mode  0  output  (fo)  by  the  output  from 
an  antenna  placed  at  the  center  of  the  array  as  sometimes  done  for  the  Watson- Watt 
approach.  If  this  were  done,  the  maximum  usable  array  size  for  the  eight-element  Butler 
matrix  approach  would  then  exceed  that  of  the  four-element  approach  as  well  as  both  the 
four  and  eight-element  Watson- Watt  approaches. 


2.3  The  Single  Signal  Maximum  Likelihood  Approach 

A  successful,  albeit  often  computationally  intensive  and  hardware  expensive,  approach 
to  estimation  is  based  on  the  maximum  likelihood  method.  Essentially  the  idea  is  to  find 
the  most  likely  state  of  a  signal  process  given  a  set  of  measurement  observations  made 
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of  this  process  plus  additive  noise.  Assuming  the  noise  has  a  Gaussian  distribution,  and 
a  measurement  of  N  sensor  outputs  is  available,  then  the  associated  probability  density 
function  is  given  by  [2] 

f(x)  =  _ I _ g-traCe((x-m)«C-i(x-m))  /.-N 

where  the  superscript  H  denotes  the  conjugate-transpose  operation,  x  is  the  noise  cor¬ 
rupted  data  vector  defined  by 

Xo 

Xl 

Xit  =  X2  (48) 

_  XN-1  _ 

m  is  the  signal  (or  noise-free)  data  vector  and  has  the  same  form  as  x,  and  C  is  the 
NxN  noise  covariance  matrix  which  describes  the  correlations  among  sensors.  The  sensor 
covariance  matrix  is  assumed  to  be  known  either  through  measurement  or  modeling. 

Under  normal  circumstances,  the  only  unknown  in  (47)  is  the  signal  data  vector  m. 
Since  the  Watson- Watt  and  Butler  matrix  approaches  make  the  assumption  that  only 
a  single  signal  is  received  at  a  particular  frequency,  the  same  assumption  is  made  here. 
Hence,  the  elements  of  the  signal  vector  m  will  be  of  the  form 

m„  =  se^¥8i“(^+^+’’)  for  n  =  0,l...,iV-l  (49) 

The  signal  vector  m  in  the  single  signal  case  is  a  function  of  the  signal  amplitude  s  and 
the  bearing  (f).  Estimating  these  two  values  can  be  accomplished  by  finding  the  choice 
of  s  and  0  which  lead  to  a  maximum  value  of  /(x).  By  inspection,  this  is  equivalent  to 
minimizing  the  expression 

=  trace((x- m)'*^C”^(x- m))  =  (x  -  m)^C“^(x  -  m)  (50) 

Defining  the  steering  vector  (or  array  response  vector)  as  e  =  m/{sy/N)  and  rewriting 
the  above  expression  in  terms  of  e  then 

=  (x-se)^C~^(x- se) 

=  x^C"^x  -  sx^C’^e  -  s*e^C~^x -I- s*se^C“^e  (51) 
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Setting  up  the  differential  equation 


yields 


e^C-^x 

eHc-ie 


which  is  the  value  of  s  which  minimizes  for  a  given  value  of  (p.  Plugging  this  value  back 
into  (51) 

(54) 


(54) 

There  is  no  obvious  way  to  directly  compute  the  value  of  0,  hence  (f>  is  estimated  by 
finding  the  value  which  minimizes  using  the  expression  shown  here,  or  equivalently  by 


maximizing 


S{cl>)  = 


e^C-^xx^C-^e 

e^C-ie 


since  the  first  term  in  (54)  is  constant  with  respect  to  (j). 


The  main  advantages  of  the  maximum  likelihood  (ML)  approach  proposed  here  is  that 
it  is  optimum  for  the  case  of  estimating  the  bearing  of  a  single  signal  in  Gaussian  noise  from 
a  single  measurement  sample,  it  is  also  unbiased  (no  need  for  a  look-up  table  for  correction 
purposes),  and  it  does  not  suffer  from  the  problem  of  the  channel  voltages  going  to  zero 
for  certain  combinations  of  signal  bearing  and  antenna  array  size.  The  disadvantages  are 
that  N  receiver  channels  are  required  for  N  antennas  and  it  is  a  narrowband  technique 
(the  steering  vector  e  is  a  function  of  the  signal  wavelength  or  frequency)  which  adds 
further  complexity  to  the  processing. 
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3.0  ERRORS  DUE  TO  NOISE  AND  INTERFERENCE 


In  this  section,  the  effects  of  internal  noise,  external  noise,  co-channel  signals  and  mul¬ 
tipath  propagation  are  considered.  These  sources  of  error  are  grouped  into  two  categories, 
namely,  internal  noise  and  external  interference.  To  simplify  the  discussion,  the  antenna 
patterns  are  assumed  to  be  ideal.  The  results  derived  based  on  these  assumptions  are 
tested  later  through  simulation. 

The  approximations  of  interest  are  taken  from  the  derivations  performed  previously 
which  make  the  assumption  that  there  is  only  one  signal  present  and  that  there  is  no 
noise.  The  reference  voltage  for  both  the  Watson- Watt  and  Butler  approaches  in  this 
case  is  then  given  by 

Vre/  ~  OfS 

where 

/Oirr  iTrn  \ 

a  =  2  cos  sin(0+ — -f  i?)j  =  2/?  (57) 

Although  a  examination  of  Figure  3  reveals  a  dependence  on  (j),  for  approximation  purposes 
it  is  simpler  to  approximate  f3  by  -I-  f3max)/2  for  the  four-element  array  and  /S^ax 
(where  Pmax  ^  Pmin)  for  the  eight-element  array.  This  leads  to  the  approximation 


where  r  <  Tmaxo- 

For  the  Watson- Watt  approach  the  relevant  channel  voltage  approximations  are 


and 


'^EW  =  jbs  sin  (f) 


(59) 


vns  =  jbs  cos  <f) 

where  the  best  choice  for  b  is  given  by 

■  N 

using  the  exact  values  of  vew  and  ujvs  described  by  equations  (15)and  (23), 


(60) 

(61) 

and  equations 
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(16)  and  (25),  respectively.  However,  this  leads  to  a  value  for  b  which  is  dependent  on  the 
bearing  <t>.  To  remove  this  dependency,  the  median  value  is  used  instead  leading  to 


3in(^)+V2sin^^j  for  =  4 

(2  +  a/2)  sin  cos  N  =  8 


where  r  <r„ 


For  the  Butler  matrix  approach  the  relevant  approximations  are 


/o  —  ^re/  — 


/i  =  cs  COS  0  +  jcs  sin 


where  the  ideal  choice  for  c  is  given  by 


mx\ 

l/ol 


using  the  exact  values  of  /i//o  given  by  (37),  (41),  and  (43)  are  used.  Since  this  again 
results  in  a  value  which  is  dependent  on  (f),  an  approximation  for  c  is  made  by  taking  the 
median  value  evaluated  for  all  values  of  ^  yielding 


sin  j  sin 


for  =  4 


,  ('Kry/2\  .  /2nr\  „  tt  ,  /27rr  7r\ 

It  j + 2  8  It  s) 

„  .  TT  .  /27rr  .  ttN 

+  2  sm  —  sm  (  — ^  sm  -  1 

8  V  A  8/ 


for  N  =  8 


where  r  <  rmax2-  The  expression  for  c  for  an  eight-element  antenna  array  is  overly  compli¬ 
cated  and,  given  the  range  of  antenna  array  sizes  of  interest,  can  be  accurately  replaced 

c  «  ^(sin(a^)+V2sin(l9^))  (67) 


The  main  advantage  of  using  the  approximations  expressed  by  (56),  (58),  (59),  (60), 
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(62),  (63),  (64),  and  (67)  is  that  the  effects  of  noise  and/or  interference  can  be  introduced 
using  superposition.  For  example,  in  the  case  of  internally  generated  noise  plus  a  single 
interfering  signal  the  reference  voltage  would  be  represented  as 


Ure/  =  as  +  aSm  +  Tl 


(68) 


where  Sm  is  the  complex  amplitude  of  the  interferer  and  n  is  the  complex  noise  value. 
Similar  modifications  would  also  be  made  to  (59),  (60),  (63),  and  (64). 


Regardless  of  the  conditions,  the  bearing  for  the  Watson- Watt  approach  is  estimated 
using 


^  =  arctan 


( 


real{-jv/^s/Vref}  J 


fxea\{-3VEwv*rtf^\ 
\  real{-jujvs  J 


(69) 


and  the  bearing  for  the  Butler  matrix  approach  is  estimated  using 


arctan 


/imag{/i//o}\ 
V  real{/i//o}  / 


arctan 


( imag{/i/o*}\ 
V  real{/i/o*}  / 


(70) 


For  the  Watson- Watt  approach,  in  the  presence  of  noise  and/or  interference  the  quantities 
VEw/vref  and  VNs/vref  are  no  longer  guaranteed  to  be  real-valued,  hence  the  reason  for  the 
real{-}  operator. 


3.1  The  Effect  of  Internal  Noise 

The  effect  of  internal  noise  is  introduced  into  (56),  (59),  (60)  (63),  and  (64)  according 
to 


^re/  — 

as +  n 

(71) 

Vew  = 

jbs  sin  4>  +  riEw 

(72) 

vns  = 

jbs  cos  (f>  -f  tins 

(73) 

fo  = 

as +  n 

(74) 

fl  = 

cs  cos  (f) +  jcs  sin  <f>  +  ni 

(75) 
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where  n  is  the  complex  noise  amplitude  in  the  reference  channel,  uew  tins  s-re  the 
complex  noise  amplitudes  in  the  directional  channels  of  the  Watson- Watt  system,  and 
rii  is  the  complex  noise  amplitude  of  noise  in  the  mode  1  channel  of  the  Butler  matrix 
system.  Since  each  receiver  channel  is  assumed  to  have  the  same  characteristics,  the  RMS 
noise  values  for  all  channels  are  assumed  to  be  the  same,  that  is 

E{|n|2}  =  E{\nEwf}  =  (76) 


3.1.1  Internal  Noise  and  the  Watson- Watt  Approach 

Given  noise  corrupted  measurements  of  the  channel  voltages,  the  estimated  signal 
bearing  using  the  Watson- Watt  approach  is  given  by 


^  =  arctan 


/ real{(&ssin 0  —  jngty)(os*  -h  n*)}' 


\real{(6scos(^  —  jnjvs)(os*  +  n*)} 

(a6|sp  sin (^  +  6 real{sn*}  sin -  imag{os*ngwr  +  ng^^n*} \ 
a6|sP  cos  (j)  -f-  6real{sn*}  cos ^  —  imag{as*njvsr  -I-  n^vsn*}  / 

Assuming  that  the  signal-to-noise  ratio  is  greater  than  a  few  dB,  and  o,  6  >  0,  then  the 
pure  noise  terms  in  the  above  expression  will  be  insignificant  and  the  expression  simplifies 
to 


-  _  f  (|sp  -I-  real{sn*/o})  sin^  —  imag{s*n£;i^/6}\ 

arctan  ^  ].eal{sn*/a})  cos^  -  imag{s*njvs/&}  / 


The  bearing  error  is  given  hy  $  —  (j>.  Taking  advantage  of  the  relationship 

tan((^  —  4>)  = 

then 

^  —  <j)  =  arctan 


(78) 


tan^  —  tan(^ 

1  -1-  tan^  tan^ 

(79) 

(  tan^  —  tan0  ^ 

Vl  +  tan^tan^/ 

(80) 
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and  incorporating  (78)  the  result  (after  removing  any  pure  noise  terms  as  was  done  pre¬ 
viously)  is  given  by 


-  ,  ^  /  imag{s*njv5/6}  sin  0  -  imag{s*n£n^/6}  cos  ^  \ 

0-0  =  arctan  - tt - ^ — .  , - — - 

\|sp  -f  real{sn*/a}  —  imag{s*njvs/o}  cos0  —  imag|s*n£w^/o}  sm0y 

(81) 

One  last  simplification  can  be  performed  by  limiting  the  results  to  bearing  errors  of  less 
than  10®  so  the  approximation  tan(0  —  0)  w  0  —  0  can  be  used,  and  also  by  recognizing 
that  for  positive  signal-to-noise  ratios  (particularly  for  levels  where  0  —  0  <  10°)  then 
|sp  »  imag{s*n/a}  cos  0,  |sp  >•  imag{s*n^’s/6}cos0  and  |sp  ^  imag{s*njsw/^}  sin  0. 
The  appropriate  simplified  error  expression  is  then  given  by 


0-0 


imag{s*n^5/6}  sin  0  -  imag{s*n£;vv/&}  cos  0 


Using  the  last  approximation  to  represent  the  bearing  error,  the  variance  of  the  error 
can  be  found  by  taken  the  expectation  of  the  squared  error  giving 


imag{(s*)^nA/’5n£’iv/6^}  sin  0  cos  0 


(imag{s*n£;vv/^>}^  cos  0)^ 


=  |s|  |imag{s*njv5/6}^}  sin^  0  -I-  |s|  ^E  |imag{s*n£;M^/6}^|  cos^  0 


where  it  has  been  assumed  that  the  signal  amplitude  is  relatively  constant,  the  noise  is 
complex  white  Gaussian  noise  so  that 

{imag{s*nivsP}  =  |sp£?  {imag{nivsp}  =  (84) 

E  {imag{s*n£;vvP}  =  |spE  {imag{n£;H^}2}  =  i|s|V  (85) 


and  where  is  the  noise  power. 
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Defining  the  signal-to-noise  ratio  as 


snr 


5| 


(86) 


then  the  RMS  bearing  error,  which  is  the  square-root  of  the  variance,  can  be  written  as 


\h\\/2snr 

Using  the  values  for  h  defined  earlier,  then  the  RMS  bearing  error  is  given  by 

1 


(87) 


sin  +  -s/Isin  \y/2snr 


\  1(2  -I-  \/2)  sin  cos  |)  -I-  v^sin  sin  |)  \y/2snr 


for  N  —  A 


for  AT  =  8 


(88) 

The  RMS  bearing  error  expressions  are  somewhat  complex,  but  for  r  A/27r  they  simplify 
to 

f  \ 

for  iV  =  4 


E{{<l>  -  <l>fY 


2\^  _ 


A'Kry/2snr 

X 

[  87rr  cos  |V2snr 


(89) 


for  AT  =  8 


Examining  the  above  expressions,  noise  related  bearing  errors  can  be  reduced  by  in¬ 
creasing  the  signal-to-noise  ratio.  Figure  12  illustrates  this  effect  for  four  and  eight-element 
arrays  showing  the  theoretical  results  using  (88)  and  the  results  based  on  computer  simu¬ 
lations  as  a  function  of  the  signal-to-noise  ratio  defined  by  (86).  For  the  simulated  results, 
1000  measurements  for  each  SNR  value  (0  to  40  dB  in  1  dB  increments)  were  used  to 
generate  the  corresponding  RMS  error  value.  Bias  errors  were  removed  before  performing 
the  RMS  calculation.  As  can  be  seen,  there  is  very  good  agreement  between  the  theoret¬ 
ical  and  simulated  results  even  for  RMS  errors  greater  than  10®.  The  eight-element  array 
also  has  approximately  half  the  RMS  error  for  a  given  SNR  than  the  four-element  array, 
as  predicted  by  (89). 

The  bearing  errors  can  also  be  reduced  by  increasing  the  size  of  the  array  up  to 
approximately  r  =  0.3A,  beyond  which  the  error  increases.  This  is  illustrated  in  Figure 
13  which  shows  the  theoretical  and  simulated  results  as  a  function  of  r  and  three  diflFerent 
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(a) 


(b) 

Figure  12:  Effect  of  the  internal  noise  level  on  the  Watson- Watt  approach  using  (a)  a 
four-element  and  (b)  an  eight-element  circular  antenna  array  with  a  radius  r  =  0.2A  and  a 
signal  bearing  of  ^  =  0°.  The  simulated  results  are  shown  by  the  solid  and  the  theoretical 
results  are  shown  by  the  dashed  line. 
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(a) 


(b) 

Figure  13:  Effect  of  the  antenna  array  radius  on  the  Watson-Watt  approach  in  the 
presence  of  internal  noise  (snr  =  10  dB)  using  (a)  a  four-element  and  (b)  an  eight-element 
circular  antenna  array.  The  theoretical  results  are  shown  by  the  dashed  line  and  the 
simulated  results  are  shown  by  the  solid  lines  for  signal  bearings  of  (i)  <f>  —  'd  =  0°,  (ii) 
—  =  22.5°,  and  (iii)  <f)  —  'd  =  45°. 
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signal  angles  for  four  and  eight-element  antenna  arrays.  For  the  simulated  results,  1000 
measurements  for  each  value  of  r  ranging  from  0.04A  to  0.65A  in  steps  of  O.OIA  were 
used  to  generate  the  corresponding  RMS  error  values.  The  sign  reversal  in  the  reference 
channel  for  r  >  Tmaxo  was  corrected  to  yield  the  bearing  in  the  correct  quadrant.  The 
simulated  results  illustrate  that  accuracy  is  also  a  function  of  the  signal  bearing.  The 
theoretical  results  provide  a  good  estimate  of  the  median  value  of  the  RMS  errors  except 
for  r  «  Tmaxo  where  the  signal  part  of  the  reference  voltage  v^f  becomes  sufficiently  small 
that  it  is  overpowered  by  noise  resulting  in  180®  ambiguity  errors  (i.e.  a  0  which 
violates  the  assumption  made  in  deriving  (78)  and  the  following  expressions).  If  a  central 
antenna  were  used  to  measure  then  the  simulated  results  would  more  closely  follow 
the  theoretical  results  in  the  region  r  «  rjnaxo- 


3.1.2  Internal  Noise  and  the  Butler  Matrix  Approach 


Given  noise  corrupted  measurements  of  the  channel  voltages,  the  estimated  signal 
bearing  using  the  Butler  matrix  approach  is  given  by 


^  =  arctan 


(■ 

( 


imag{ (cs  cos  (f)  -I-  jcs  sin  (j)  -f-  ni) (as*  +  n*)}  A 
real{ (cs  cos  <t>  +  jcs  sin  <f)-¥n\) (as*  +  n*) }  j 

(acjsp  -I-  real{csn*})  sin  <f>  +  imag{as*ni  -i-  csn*  cos  (j)  -1-  nin*} 


=  arctan  |  y  i-,  - 1, —  j,  —  t  ■  - o<. —  -x  ■  —  — -r  ■  -i--  j  /hQS 

(ac|sP  +  real{csn*})  cos  0  -I-  real{as*ni  -f  nin*}  —  imag{csn*sin(^}/ 

Assuming  that  the  signal-to-noise  ratio  is  greater  than  a  few  dB,  then  the  pure  noise 
terms  in  the  above  expression  will  be  insignificant  and  the  expression  simplifies  to 


<!> 


^  ,  (acs^  +  real{csn  })sm(i  +  imag{as  ni-fcsn  cos^)  \  ,  . 

arctan  I  7 — ,  — ,,  (91) 

(ac|s|'^  +  real{csn*})  cos  (j)  -f  real{as*ni}  —  imag{csn*  sin  J 


( 


Taking  advantage  of  the  relationship  expressed  in  (80),  the  bearing  error  can  be  ex¬ 
pressed  as 


6  —  6  =  arctan  I  — oy-  -  /  ~j  - l-  --iz-j — r-  •  .  .  . 


imag{ sn*/ a}  —  real{s*ni/c}  sin  ^  -I-  imag{s*ni/c}  cos  (j)  ' 

-I-  real{sn*/a}  +  real{s*ni/c}  cos  (f>  -I-  imag{s*ni/c}  sin  <f>^ 


where  pure  noise  terms  have  again  been  removed. 

As  was  done  for  the  Watson- Watt  approach,  one  last  simplification  can  be  performed 
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by  limiting  the  results  to  bearing  errors  of  less  than  10°  so  the  approximation  tan(^  — 
can  be  used,  and  also  by  recognizing  that  for  positive  signal-to-noise  ratios 
(particularly  for  levels  where  ^  —  10°)  then  |sp  '>  |s*n/ap  and  |sp  '>  js*ni/cp.  The 

appropriate  simplified  error  expression  is  then  given  by 

-  ,  imag{sn*/a}  —  real{s*ni/c}sin^  +  imag{s*ni/c}cos^ 

^-■i>  =  - 

Using  the  above  approximation  to  represent  the  bearing  error,  the  variance  of  the  error 
can  be  found  by  taken  the  expectation  of  the  squared  error  giving 


Eui'-m  =  e\ 

fimag{sn7oPl  ,  E,/real{ 

l 

s*ni/c}^sin^ 

/ 

,  ^  fimag{s*ni/c}2cos2 

1  1 

-  2\sP'-g?'^  c‘’ 

(94) 

where  it  has  been  assumed  that  the  signal  amplitude  is  deterministic  and  the 
complex  white  Gaussian  in  nature  so  that 

noise  is 

E  |imag{sn*}^j 

[  =  ls|2£;{imag{n*p} 

(95) 

E  |real{s*ni}^j 

[  =  |sp£'  {imag{ni}^} 

= 

(96) 

j&{imag{s*ni}^j 

[  =  |sp£'{imag{nip} 

=  jWV 

(97) 

and  where  is  the  noise  power. 

Using  the  signal-to-noise  ratio  defined  previously  in  (86),  the  RMS  bearing  error,  which 
is  the  square  root  of  the  variance,  can  be  written  as 
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(a) 


(b) 

Figure  14:  Effect  of  the  internal  noise  level  on  the  Butler  matrix  approach  using  (a)  a 
four-element  and  (b)  an  eight-element  circular  antenna  array  with  a  radius  r  =  0.2A  and  a 
signal  bearing  of  0  =  0®.  The  simulated  results  are  shown  by  the  solid  and  the  theoretical 
results  are  shown  by  the  dashed  line. 
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Applying  the  definition  for  a  given  by  (58)  and  the  definition  for  c  given  by  (67)  the  RMS 
bearing  error  expression  can  be  expanded  to  become 


For  small  array  sizes  where  r  A/27r  then  a  «  iV  and  c  w  Nttt/X  and  the  RMS  bearing 
error  expression  can  also  be  written  in  a  simpler  form,  namely, 


1  / 

N^/2anr\  irvj 


(100) 


Examining  the  RMS  bearing  error  expression,  noise  related  bearing  errors  can  be 
reduced  by  increasing  the  signal-to-noise  ratio.  Figure  14  illustrates  this  effect  for  an 
eight-element  array  with  r  =  0.2A  using  the  same  data  used  to  generate  Figure  12  for  the 
Watson-Watt  approach.  Bias  errors  were  removed  before  performing  the  RMS  calculation. 
As  can  be  seem,  there  is  very  good  agreement  between  the  theoretical  and  simulated  results 
even  for  RMS  errors  greater  than  10°. 

Although  less  obvious,  the  bearing  errors  can  also  be  reduced  by  increasing  the  size 
of  the  array  up  to  r  <  0.19A,  beyond  which  the  error  increases  substantially  as  noise 
begins  to  dominate  the  phase  of  the  mode  0  output  /q.  Since  the  bearing  is  derived  from 
the  phase  of  fi/fo,  these  phase  errors  translate  directly  into  bearing  errors.  Figure  15 
illustrates  this  effect  for  an  eight-element  array  for  r  «  rmaxo  using  the  same  data  that 
was  used  to  generate  Figure  13  for  the  Watson-Watt  approach.  The  RMS  errors  were 
found  to  be  relatively  independent  of  direction  so  only  one  signal  bearing  was  simulated. 
In  theory  the  large  errors  for  r  «  rmaxo  could  have  been  eliminated  by  replacing  the  mode 
0  output  fo  by  the  output  from  a  central  antenna.  The  large  errors  as  r  — >  rmax2  would  be 
unaffected  by  this  modification  and  would  limit  the  maximum  usable  size  of  the  antenna 
array.  The  agreement  between  the  theoretical  and  simulated  results  is  very  good  in  this 
case. 
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(a) 


(b) 

Figure  15:  Effect  of  the  antenna  array  radius  on  the  Butler  matrix  approach  in  the 
presence  of  internal  noise  {snr  =  10  dB)  using  (a)  a  four-element  circular  antenna  array 
with  signal  bearings  of  (i)  =  0°,  (ii)  ^  =  22.5®,  and  (iii)  <j)-'d  =  45°,  and  (b)  an 

eight-element  circular  antenna  array  with  a  signal  bearing  ^  =  0®  (results  for  other 

signal  bearings  are  the  same).  The  theoretical  results  are  shown  by  the  dashed  line  and 
the  simulated  results  are  shown  by  the  solid  lines. 
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3.1.3  Internal  Noise  and  Comparisons  with  Maximum  Likelihood  Perfor¬ 
mance 

The  performance  of  the  Watson- Watt  (worst  case),  Butler  matrix,  and  ML  approaches 
is  compared  in  Figures  16  and  17  using  the  same  data  used  to  generate  Figures  12  and  13 
respectively.  Figure  16  shows  the  Watson- Watt  and  Butler  matrix  results  are  similar  and 
significantly  better  than  the  Maximum  Likelihood  approach  when  using  an  array  size  of 
r  =  0.2A.  This  somewhat  surprising  result  is  due  to  the  fact  that  the  Maximum  Likelihood 
approach  uses  an  eight-channel  receiver  with  noise  introduced  into  all  eight  channels.  The 
Watson- Watt  and  Butler  matrix  approaches  use  a  lesser  number  of  receivers  (three  and 
two-channel  receivers,  respectively)  resulting  in  a  correspondingly  lesser  amount  of  noise 
being  introduced.  This,  in  turn,  results  in  improved  performance. 

Figure  17  highlights  the  differences  between  all  three  approaches  particularly  as  the 
radius  of  the  array  approaches  rmaxo  =  0.35A  for  the  four-element  array  and  rmaxo  =  0.38A 
for  the  eight-element  array.  At  this  radius  the  gain  of  the  reference  channel  goes  to 
zero  and  the  bearing  estimates  become  extremely  susceptible  to  noise  errors.  In  the  case 
of  the  Watson- Watt  approach,  the  resultant  phase  errors  only  become  large  enough  at 
r  =  rmaxo  ~  0.03A  to  cause  180“  ambiguity  errors  which,  in  turn,  dramatically  increases 
the  RMS  bearing  error  level  (note  that  for  lower  SNR  ambiguity  errors  occur  at  a  smaller 
value  of  r  while  for  higher  SNR  ambiguity  errors  occur  at  a  larger  radius  but  not  exceeding 
rmaxo)-  For  the  Butler  matrix  approach,  the  effect  of  phase  errors  in  the  reference  channel 
is  more  gradual  and  begins  at  a  smaller  value  of  r  so  that  performance  is  always  worse 
than  that  of  the  Watson- Watt  approach.  The  ML  approach  does  not  suffer  zero  gain 
problems,  but  for  the  four-element  array,  ambiguities  in  the  bearing  estimates  occur  for 
r  =  rmaxo  resulting  in  large  bearing  errors.  For  the  eight-element  array,  there  are  no  such 
problems  so  that  its  performance  continuously  improves  as  a  function  of  the  array  radius 
up  to  r  =  0.7A  (not  shown),  and  it  also  attains  the  lowest  RMS  error  for  r  >  0.3A  despite 
the  extra  noise. 

The  problem  due  to  the  reference  channel  could  be  overcome  by  using  a  dedicated 
antenna  positioned  at  the  center  of  the  antenna  array  as  mentioned  previously.  In  the  case 
of  the  Watson- Watt  approach,  performance  would  then  be  limited  by  the  gain  becoming 
too  low  in  the  two  directional  channels  (i.e.,  vew>  '^^ns  — >•  0)  which  occurs  when  r  rmaxi 
(0.5A  for  the  four-element  array  and  0.54A  for  the  eight-element  array).  In  the  case  of  the 
Butler  matrix  approach,  performance  would  be  limited  by  the  gain  becoming  too  low  in 
the  mode  1  channel  (i.e.,  /i  -¥  0)  which  occurs  when  r  -¥  Tmax^  (0.5A  for  the  four-element 
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(a) 


(b) 

Figure  16:  Effect  of  internal  noise  on  the  Watson- Watt,  Butler  matrix,  and  ML  ap¬ 
proaches  using  (a)  a  four-element  and  (b)  an  eight-element  circular  array  with  a  radius 
r  =  0,2A  and  a  signal  bearing  of  ^  =  0°. 
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Figure  17:  Effect  of  the  antenna  array  radius  on  the  performance  of  the  Watson- Watt, 
Butler  matrix,  and  ML  approaches  in  the  presence  of  internal  noise  {snr  =  10  dB)  using  (a) 
a  four-element  circular  array  with  a  signal  bearing  of  i?  =  45®  and  (b)  an  eight-element 
circular  array  with  a  signal  bearing  of  ^  =  0®. 
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array  and  0.6A  for  the  eight-element  array). 


3.2  The  Effect  of  External  Interference 

As  mentioned  previously,  the  term  “external  interference”  is  used  to  include  external 
noise,  multipath,  and  co-channel  interference.  Although  these  are  different  phenomena, 
they  can  be  handled  in  a  very  similar  way  since  they  can  all  be  modeled  as  a  number 
of  interfering  signals.  The  differences  in  these  phenomena  are  taken  care  of  by  adjusting 
the  numbers  of  interfering  signals,  the  statistical  distribution  of  these  signals,  and  their 
temporal  and  spatial  correlations. 

For  this  study,  external  noise  was  modeled  as  a  very  large  number  of  interfering  signals 
arriving  from  a  wide  range  of  directions.  The  advantage  of  this  approach  is  that  it 
automatically  takes  into  account  the  correlations  between  the  noise  measured  at  different 
sensors  -  correlations  which  can  be  significant  when  the  spacings  are  less  than  A/2.  The 
complex  amplitudes  of  these  signals  were  assumed  to  have  a  random  Gaussian  distribution 
and  to  be  uncorrelated  from  measurement  to  measurement  (temporal  decorrelation),  and 
also  to  be  uncorrelated  if  the  DF  system  was  moved  from  site  to  site  (spatial  decorrelated). 
The  bearings  of  the  interfering  signals  were  assumed  to  be  uniformly  distributed  360°  in 
azimuth  and  also  temporally  and  spatially  uncorrelated.  Note  that  this  best  models  man¬ 
made  noise  which  has  terrestrial  origins,  but  not  isotropic  or  cosmic  noise  which  has  an 
elevation  bearing  component  as  well.  The  modification  for  isotropic  noise  is  discussed 
later  on. 

Co-channel  interference  was  also  modeled  in  the  same  way  as  noise  except  the  numbers 
of  interfering  signals  were  assumed  to  be  small  (i.e.  one  to  several  signals). 

For  modeling  multipath  signals,  there  were  several  key  differences.  The  numbers  of 
interfering  signals  were  assumed  to  range  from  a  small  value  for  specular  multipath  to  a 
large  value  for  diffuse  multipath.  The  complex  amplitudes  of  these  signals  were  assumed 
to  have  a  random  Gaussian  distribution  when  measured  from  site  to  site  (spatially  un¬ 
correlated),  but  to  be  invariant  with  time  when  measured  for  at  a  given  site  (temporally 
correlated).  The  bearings  were  assumed  to  have  distributions  which  ranged  from  uniform, 
for  sources  of  multipath  close  to  the  receiver,  to  concentrated  in  the  direction  of  <f>,  for 
sources  of  multipath  far  from  the  receiver  (and  close  to  the  transmitter).  Other  distribu¬ 
tions  are  possible  but  were  not  investigated  in  this  study.  The  signal  bearings  were  also 
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assumed  to  be  spatially  uncorrelated  but  temporally  correlated.  The  assumption  of  tem¬ 
poral  correlation  will  not  necessarily  be  true  for  all  cases  of  multipath  (e.g.  HF  skywave 
signals  where  path  length  differences  can  be  large  enough  that  decorrelation  does  occur), 
however  the  following  results  are  not  changed  as  long  as  there  is  either  spatial  or  temporal 
decorrelation. 

For  any  of  these  three  types  of  interference,  the  effects  of  the  interfering  signals  can 
be  introduced  into  (56),  (59),  (60),  (63),  and  (64)  by  adding  the  new  signals  according  to 


Vref  = 

aJ2si 

i=0 

(101) 

Vew  = 

jb'^  Si  sin (t)i 

1=0 

(102) 

vns  = 

jh'^Si  cos  (j)i 

i=0 

(103) 

fo  = 

1=0 

(104) 

h  = 

c  ^  Si  cos  (j>i  +  j  Si  sin  (j)i 

1=0 

(105) 

where  the  summation  limit  for  i  is  not  shown  but  assumed  to  be  chosen  appropriate  to 
the  interference  mechanism,  so  =  s  is  the  complex  amplitude  of  the  signal  of  interest  and 
^0  =  ^  is  corresponding  bearing  (the  subscript  0  was  added  for  convenience),  and  si,  S2,  — 
are  the  complex  amplitudes  of  the  interfering  signals  and  are  the  corresponding 

bearings. 


3.2.1  Interference  and  the  Watson- Watt  Approach 


Given  measurements  of  the  channel  voltages  corrupted  by  external  interfering  signals, 
the  estimated  signal  bearing  using  the  Watson- Watt  approach  is  given  by 


<j>  —  arctan 


/ real{(5  Ei=o  Sj  sin  (f>i)  (a  Ei=o  <)}  \ 
Vreal{(6  Ei=o  Si  cos  <^i)(a  Ei=o  «* )}  / 


=  arctan 


/  Ei=o  Ek=o  rea.l{sisl}  sin  (l>i  \ 
VEi=o  Efc=o  real{sis;^}  cos  4)i  J 


(106) 
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where  the  assumption  that  a,b  >  0  has  been  made.  Using  the  relationship  expressed  in 
(80),  and  after  some  rearranging  the  bearing  error  can  be  expressed  as 


f  Ei=o  Ek=o  real{sfs;^}  sin(^i  -  ^p)  N 
VEi=o  Efc=o  real{s<s;^}  cos{<f>i  -  <f>o)  J 


(107) 


Limiting  the  results  to  bearing  errors  of  less  than  10°,  then  the  approximation  w 

^  —  <i>  will  be  sufficiently  accurate  for  the  purposes  here.  Additionally,  in  order  for  the  10° 
limit  to  be  true  then  |so|  ^  |st|  for  i  =  1, 2, ...  which  leads  to  the  simplified  expression 


Ei=i  real{sisS}  sin(^i  -  ^o) 


(108) 


For  a  specific  measurement,  given  that  the  signal  parameters  are  all  known,  then  (108) 
can  be  used  to  determine  the  expected  bearing  error.  This  can  be  useful  when  exploring 
specific  effects  of  co-channel  interference  or  multipath  through  computer  simulation.  In 
most  real-life  applications,  however,  the  signal  parameters  will  not  be  known  so  it  is  of 
more  interest  to  statistically  quantify  the  errors.  The  most  useful  quantity  in  this  regards 
is  the  RMS  value.  In  the  following  derivation  of  the  RMS  value,  the  statistical  results  are 
based  on  measurements  taken  at  different  time  instances  and  different  measurement  sites 
(while  keeping  the  parameters  of  the  signal  of  interest  fixed),  hence  the  temporal  and/or 
spatial  decorrelation  properties  of  the  interfering  signal  parameters  can  be  utilized. 

The  RMS  value  can  be  determined  by  first  finding  the  variance  of  the  error  which  is 
given  by 


E{{^  -  0)2} 


p  f  (Ei=i  realfsjsg}  sin(0i  -  0o))2  } 

I  kor  j 

E  {E,=i  real{sisS}2  s\r^{(j)i  -  0o)} 

SF 

E  {Ei=i  rea\{sisl}^}  E  {sin2(0,  -  0o)} 


4|s|2 


(109) 


where  s  =  sq,  is  the  total  interference  signal  power  given  by  -i- 


39 


is  the  interference  distribution  factor  which  will  be  discussed  shortly,  and  the  relationships 


and 


E{sin\(l>i-<i>o)} 


J3^{real{5is;}2}  =  ^|s|V 

i=l  ^ 


(110) 


sin^{^j  -  ^o)#i  +  (1-Q^)£{sin^(^-^)}  (111) 


(112) 


were  used.  The  factor  depends  on  the  angular  distribution  of  the  interference  source 
and  for  a  uniform  distribution  over  360°  (represented  by  the  first  term  in  (111))  then 
0^  =  1,  and  for  a  tight  distribution  centered  around  0  (represented  by  the  second  term  in 
(111))  then  «  0.  For  external  noise  and  co-channel  interference  =  1  can  normally 
be  assumed.  For  multipath,  is  interpreted  as  the  power  ratio  between  the  part  of 
the  multipath  generated  as  a  uniform  distribution  and  the  total  multipath  power.  The 
value  is  dependent  on  the  signal  environment.  For  example,  given  a  large  number  of 
small  (relative  to  the  wavelength)  scattering  sources  of  multipath  and  a  terrestrial  signal 
path,  most  of  the  multipath  is  generated  near  the  receiver  (uniform  distribution),  near 
the  transmitter  and  along  the  signal  path  (tight  distribution)  leading  to  a  value  «  0.5 
[3].  In  another  example,  given  a  large  number  of  large  reflecting  sources,  the  multipath 
region  is  much  larger  and  more  circular  than  the  scattering  case  leading  to  a  value  «  1 
[3].  One  final  example  is  a  calibrated  airborne  surveillance  system,  where  most  of  the 
multipath  will  be  produced  near  the  transmitter  resulting  in  a  value  of  «  0. 


Finally,  defining  the  signal-to-interference  power  ratio  as 


(113) 


(where  sir  is  replaced  by  snr  for  external  noise  calculations)  then  the  RMS  bearing  error 
is  given  by 


(114) 


An  interesting  feature  of  this  expression  is  that  unlike  the  expressions  (88)  and  (99), 
developed  to  quantify  the  effects  of  internal  noise,  the  size  of  the  array  and  number  of 
antennas  are  not  factors,  only  the  signal-to-interference  power  ratio.  Additionally,  these 
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(a) 


(b) 

Figure  18:  Effect  of  external  terrestrial  interference  on  a  Watson- Watt  system  using  (a) 
a  four-element  and  (b)  an  eight-element  circular  antenna  array  with  a  radius  r  —  0.2A  and 
a  signal  bearing  =  The  simulated  results  are  shown  by  the  solid  and  the  theoretical 
results  are  shown  by  the  dashed  line. 
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(a) 


(b) 

Figure  19:  Effect  of  the  antenna  array  radius  on  the  Watson- Watt  approach  in  the 
presence  of  external  terrestrial  interference  (sir  =  10  dB)  using  (a)  a  four-element  circular 
antenna  array  with  signal  bearings  of  (i)  <j)—'d  —  0°,  (ii)  t?  =  22.5°,  and  (iii)  =  45°, 
and  (b)  an  eight-element  circular  antenna  array  with  a  signal  bearing  (^  —  =  0°  (results 

for  other  signal  bearings  are  the  same).  The  theoretical  results  are  shown  by  the  dashed 
line  and  the  simulated  results  are  shown  by  the  solid  lines. 
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expressions  suggest  that  the  low  gain  problem  encountered  in  the  presence  of  internal 
noise  when  r  «  Tmaso  ot  r  fa  Tmaxi  will  not  be  a  factor  since  the  gain  of  the  signal  and 
interference  will  both  be  effected  in  the  same  way  leaving  the  sir  unchanged  (assuming  of 
course  that  internal  noise  is  absent  -  an  unrealistic  assumption  in  practical  applications 
but  useful  for  identifying  the  causes  of  error  as  done  here).  To  test  these  conclusions, 
simulations  were  carried  out  in  the  same  way  as  was  done  for  the  results  displayed  in 
Figures  12  and  13)  except  the  internal  noise  was  replace  by  external  terrestrial  noise.  The 
results  are  shown  in  Figures  18  and  19. 

The  simulated  results  indicate  that,  as  predicted,  the  RMS  errors  are  relatively  inde¬ 
pendent  of  the  size  of  the  antenna  array  except  for  radii  approaching  the  critical  values 
I"  —  Tmaxo  and  r  =  Tmaxi  •  The  explanation  for  the  failure  of  the  theoretical  expressions  to 
explain  these  critical  array  radii  is  due  to  the  incorrect  assumption  that  the  appropriate 
channel  gains  go  to  zero  at  r  =  Tmaxo  and  r  =  rman  for  all  signal  directions.  This  is  readily 
apparent,  for  example,  when  one  considers  that  cx:  ^  and  observes  the  variations  in 
P  shown  as  a  function  of  signal  direction  in  Figure  3  for  a  four-element  antenna  array. 
The  result  is  that  the  gain  for  the  interference  arriving  from  bearings  different  from  the 
signal  will  not  necessarily  go  to  zero  when  the  signal  gain  does  so  that  the  sir  decreases 
significantly  and  large  errors  result.  The  same  problem  occurs  for  r  =  rmoxi  for  both  the 
four  and  eight-element  antenna  arrays.  An  exception  is  for  r  =  rmcao  for  the  eight-element 
array  since  at  this  radius  the  gain  of  the  reference  channel  is  virtually  independent  of 
signal  direction  (see  Figure  4)  so  that  the  simulated  results  are  much  closer  to  the  theo¬ 
retical  expectations.  However,  even  in  this  special  case,  under  more  realistic  conditions, 
the  effects  of  internal  noise  would  certainly  dominate. 

The  effect  of  isotropic  noise  can  be  easily  introduced  in  the  previous  derivation  by 
replacing  every  term  sin  (f>i  by  sin  0*  cos  tpi,  and  every  term  cos  by  cos  (pi  cos  ipi  where  ipi  is 
the  elevation  angle  and  where  the  receive  antennas  are  assumed  to  have  an  omnidirectional 
gain  pattern.  The  signal  of  interest  is  still  assumed  to  have  no  elevation  so  tpo  =  0“.  For 
isotropic  noise  the  interfering  signals  will  be  uniformly  distributed  over  the  half  sphere 
formed  by  0°  <  ^  <  360°  and  0°  <  ^  <  V’/2,  hence  (112)  becomes 

E  {sin2(^i  -  (po)  cos2(^<)}  = 


/•27r  /*f  « 

/  /  sin^(^i  -  <po)  cos^('(pi)dipid(pi 

Jo  Jo 


J_  /■! 

27r 


1 

3 


(115) 
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where  =  1  was  used,  and  (114)  becomes 

EU}  - 


1 

\/6snr 


(116) 


where  snr  was  used  instead  of  sir  since  noise  is  being  discussed.  Interestingly  for  the  same 
SNR,  isotropic  noise  causes  less  error  than  terrestrial  noise  implying  that  interference  from 
high  elevation  angles  has  less  effect  than  interference  at  low  elevation  angles.  Figures  20 
and  21  show  the  same  noise  effects  as  Figures  18  and  19  except  isotropic  noise  was  used 
instead  and  only  the  results  for  the  eight-element  antenna  array  have  been  shown.  In 
comparing  the  results  of  isotropic  noise  to  terrestrial  noise,  the  same  comments  about 
the  terrestrial  noise  apply  here  except  that  the  zero  gain  problem  has  been  exacerbated. 
In  this  case,  for  high  elevation  angle  noise  interference  the  gain  does  not  go  to  zero  at 
f  =  fmaxo  ^  =  ^maxi  and  Consequently  a  greater  amount  of  interference  is  available  to 
overpower  the  signal  causing  large  errors  for  these  critical  antenna  array  sizes. 


Figure  20:  Effect  of  external  isotropic  noise  on  a  Watson- Watt  system  using  an  eight- 
element  circular  antenna  array  with  a  radius  r  =  0.2A.  The  simulated  results  are  shown 
by  the  solid  and  the  theoretical  results  are  shown  by  the  dashed  line. 
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Figure  21:  Effect  of  the  antenna  array  radius  on  a  Watson- Watt  system  using  an  eight- 
element  circular  antenna  array  in  the  presence  of  external  isotropic  noise  {snr  =  10  dB). 
The  simulated  results  are  shown  by  the  solid  and  the  theoretical  results  are  shown  by  the 
dashed  line. 


3.2.2  Interference  and  the  Butler  Matrix  Approach 


Given  channel  voltage  measurements  corrupted  by  external  interference,  the  estimated 
signal  bearing  using  the  Butler  matrix  approach  is  given  by 


arctan 


arctan 


/  imag{(c  E,=o  Sj  cos  -I-  jsj  sin  <j)i){a  Y.k=o  4)}  \ 
V  real{ (c  E<=o  Sf  cos  4>i  jsi  sin  (pi)  (a  Zk=o  4)}  / 

(Ei=o  Efc=o  real{sis]fc}  sin  (pj  +  imag{siSfc}  cos  cpj  > 
Ei=o  Efc=o  real{si5]^}  cos  -  imag{si5]J}  sin  (pi  j 


where  the  assumption  that  a, c>  0  has  been  made. 


(117) 


Using  the  relationship  expressed  in  (80)  and  then  rearranging,  the  bearing  error  can 
be  expressed  as 


(p  —  ^  =  arctan 


{ Et=o  Efc=o  real{sis;^}  sin(<^,-  -  (pg)  +  imagfsjs;^}  cos{(pi  -  0o)  \ 
VEj=o  Eifc=o  realisfsjfe}  cos{(pi  -  <^o)  -  imag{sis;t}  sin(^i  -  <^o)/ 


As  has  been  done  previously,  by  limiting  the  results  to  bearing  errors  of  less  than  10°  then 
the  approximation  tan(^  —  (p)  fa  ^  —  (p  can  be  used.  Additionally,  in  order  for  the  10° 
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limit  to  be  true  then  |so|  >  |si|  for  i  =  1, 2, which  leads  to  the  simplified  expression 


7  ,  Ei=i  real{sisS}  sm{<f)i  -  ^o)  +  imag{ SiSj} (cos ((^i  -  (po)  -  1) 

M 


(119) 


Using  the  same  line-of-reasoning  to  develop  an  RMS  bearing  error  expression  in  the 
previous  section  for  the  Watson-Watt  approach,  the  RMS  error  can  be  found  by  first 
determining  the  variance  according  to 


=  E 


(Ei=i  real{sjSo}  -  (j)o)  +  imag{siSg}(cos(^i  -  0o)  -  1))^ 


E  {St=i  real{stSo}^  sin^(<^t  -  (^p)  +  imag{stSo}^(cos(0i  -  <f>o)  -  1)^} 


(120) 


where  the  relationships  in  (110)  and  (112)  were  used,  as  well  as 


53£^{imag{sisS}2}  = 


(121) 


E  {(cos(^i  -  <^o)  -  1)^}  =  ^  ^  (cos(^i  -  (f>o)  -  l)^#i 

.+  (1  -  a2)£;{(cos(<^o  -  <^o)  -  1)')}  (122) 


(123) 


and  where  is  the  angular  distribution  factor  for  multipath  interference  discussed  in  the 
previous  section.  For  noise  and  co-channel  interference  0^  =  1. 

Finally,  using  the  definition  for  the  signal-to-interference  power  ratio  given  earlier  in 
(113),  the  RMS  bearing  error  is  given  by 


(124) 


As  in  the  case  of  (114),  the  expression  for  Watson-Watt  approach,  external  interference 
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(noise,  co-channel  interference,  or  multipath)  induced  bearing  errors  are  reduced  by  in¬ 
creasing  the  signal-to-interference  power  ratio,  but  are  unaffected  by  the  size  of  the  array. 
The  RMS  bearing  error  for  the  Butler  matrix  approach,  however,  is  twice  that  for  the 
Watson- Watt  approach.  This  is  a  consequence  of  the  fact  that  accurate  quadrature  sepa¬ 
ration  of  the  cos  (f)  and  sin  (j)  spatial  signal  components  in  the  /i  channel  is  dependent  on 
the  phase  accuracy  of  the  reference  (/o)  channel.  Any  phase  errors  in  /o  translate  directly 
into  bearing  errors.  In  the  Watson- Watt  approach  the  reference  channel  is  only  used  to 
resolve  180®  bearing  ambiguities  and  bearing  errors  are  only  introduced  when  gross  phase 
errors  in  the  reference  channel  occur. 

A  comparison  between  the  theoretical  and  simulated  results  is  shown  in  Figures  22 
and  23.  The  data  was  produced  in  the  same  manner  as  the  data  used  for  Figures  18  and 
19.  Since  the  qualitative  similarities  and  differences  between  the  simulated  and  predicted 
results  are  nearly  identical  to  those  for  the  Watson- Watt  approach  in  Figures  18  and  19, 
see  the  discussion  in  the  previous  section. 

The  effect  of  isotropic  noise  can  be  introduced  into  the  previous  derivation  by  replacing 
every  term  sin  <l>i  by  sin^jCosV’i,  and  every  term  cos<f>i  by  cos^jcos^'i  where  is  the 
elevation  angle  and  where  the  receive  antennas  are  assumed  to  be  omnidirectional.  The 
signal  of  interest  is  still  assumed  to  have  no  elevation  so  =  0°.  For  isotropic  noise  the 
interfering  signals  will  be  uniformly  distributed  over  the  half  sphere  formed  by  0°  <  ^  < 
360®  and  0®  <  ■0  <  'ip/2,  hence  (112)  becomes 

E  {(cos(0i  -  ^o)  cos(V’i)  -  1)^}  =  ^  ^  (cos(^i  -  ^o)  cos(^,)  -  1)^  cos{'ipi)d'ipid(pi 

=  I  (125) 

where  =  1  was  used.  Using  the  relationship  in  (115)  as  well,  then  (124)  becomes 

(126) 

where  the  term  sir  was  replaced  by  snr  since  noise  is  being  discussed.  Inspecting  (124) 
and  (126),  the  change  in  noise  distribution  from  zero-elevation  to  isotropic  results  in  lower 
error,  although  the  ratio  of  improvement  is  slightly  less  than  for  the  Watson- Watt  system. 
A  comparison  of  the  theoretical  and  simulated  results,  for  the  eight-element  antenna  array 
only,  is  provided  in  Figures  24  and  25  which  are  a  repeat  of  Figures  22  and  23  except 
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(a) 


(b) 

Figure  22:  Effect  of  external  terrestrial  interference  on  a  Butler  matrix  system  using  (a) 
a  four-element  and  (b)  an  eight-element  circular  antenna  array  with  a  radius  r  =  0.2A  and 
a  signal  bearing  (j)  =  0°.  The  simulated  results  are  shown  by  the  solid  and  the  theoretical 
results  are  shown  by  the  dashed  line. 
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(a) 


(b) 


Figure  23:  Effect  of  the  antenna  array  radius  on  the  Butler  matrix  approach  in  the 
presence  of  external  terrestrial  interference  {sir  =  10  dB)  using  (a)  a  four-element  circular 
antenna  array  with  signal  bearings  of  (i)  i?  =  0®,  (ii)  <f>—'d  =  22.5°,  and  (iii)  (^— i?  =  45°, 
and  (b)  an  eight-element  circular  antenna  array  with  a  signal  bearing  ^  —  t?  =  0°  (results 
for  other  signal  bearings  are  the  same).  The  theoretical  results  are  shown  by  the  dashed 
line  and  the  simulated  results  are  shown  by  the  solid  lines. 
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isotropic  noise  was  used  instead. 


In  comparing  the  results  for  isotropic  noise  to  terrestrial  noise,  the  same  comments 
about  the  terrestrial  noise  apply  here  except  that  the  zero  gain  problem,  like  the  Watson- 
Watt  approach,  has  been  exacerbated.  The  gain  for  high  elevation  angle  noise  interference  ^ 

does  not  go  to  zero  at  r  =  Tmaxo  and  r  =  r^aii  and  consequently  a  greater  amount  of 
interference  is  available  to  overpower  the  signal  causing  large  errors  for  these  critical  > 

antenna  array  sizes. 


Figure  24:  Effect  of  external  isotropic  noise  on  a  Butler  matrix  system  using  an  eight- 
element  circular  antenna  array  with  a  radius  r  =  0.2A.  The  simulated  results  are  shown 
by  the  solid  and  the  theoretical  results  are  shown  by  the  dashed  line. 
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Figure  25:  Effect  of  the  antenna  array  radius  on  a  Butler  matrix  system  using  an  eight- 
element  circular  antenna  array  in  the  presence  of  external  isotropic  noise  {snr  =  10  dB). 
The  simulated  results  are  shown  by  the  solid  and  the  theoretical  results  are  shown  by  the 
dashed  line. 

3.2.3  Interference  and  Comparisons  with  Maximum  Likelihood  Performance 

The  performance  of  the  Watson- Watt,  Butler  matrix,  and  Maximum  Likelihood  ap¬ 
proaches  is  compared  in  Figures  26-29  using  the  same  data  used  to  generate  Figures  18-21 
respectively. 

Figure  26  shows  performance  (i.e.  RMS  bearing  errors)  as  a  function  of  the  signal- 
to-interference  ratio,  and  unlike  the  internal  noise  case,  the  ML  approach  significantly 
outperforms  both  the  Watson- Watt  and  Butler  matrix  approaches.  Increasing  the  number 
of  antenna  elements  from  four  to  eight  also  improves  ML  performance  but  has  no  effect 
on  either  the  Watson- Watt  or  Butler  matrix  approaches.  The  performance  of  the  Butler 
matrix  suffers  the  most  due  to  its  increased  sensitivity  to  phase  errors  in  the  reference 
channel. 

Figure  27  illustrates  that  the  performance  of  the  Watson- Watt  and  Butler  matrix  ap¬ 
proaches  are  independent  of  the  number  of  antenna  elements  and  array  size  for  smaller 
arrays.  Performance,  however,  dramatically  worsens  (i.e.  the  RMS  bearing  errors  in¬ 
crease)  as  the  array  radius  approaches  the  critical  values  of  r^oao,  Tmaxi  (Watson- Watt) 
or  rmax2  (Butler  matrix).  Outside  these  critical  values,  the  Butler  matrix  approach  is 
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approximately  a  factor  of  two  times  worse  than  the  Watson- Watt  approach.  By  way  of 
comparison,  the  performance  of  the  ML  approach  is  considerably  better,  but  it  too  suffers 
the  same  problems  as  the  Watson- Watt  and  Butler  matrix  approaches  when  the  radius 
approaches  the  critical  values  of  rmaxo^  fmaxi  and  for  the  four-element  array.  There 
is  no  such  critical  radius  problem  for  the  eight-element  array,  however,  and  the  ML  ap¬ 
proach  works  well  even  for  array  radii  exceeding  rmax2-  Overall  the  performance  of  the 
ML  approach  was  a  factor  of  two  to  five  times  better  than  the  Watson- Watt  approach,  a 
performance  improvement  which  can  be  largely  attributed  to  the  ML  approach’s  ability 
to  take  advantage  of  the  second  order  statistical  nature  of  the  noise  (which  is  assumed  to 
be  known  or  measureable)  as  well  as  the  advantage  of  using  N  receiver  channels  instead 
of  two  or  three. 

Figures  28  and  29  illustrate  the  same  effects  as  Figures  26  and  27  except  that  isotropic 
noise  was  used  instead  and  only  the  results  for  an  eight  element  antenna  array  are  shown. 
Relative  performances  remained  unchanged  and  absolute  performance  was  improved  for 
small  array  sizes,  but  worse  for  an  array  radius  approaching  rmaxo-  The  poor  performance 
for  array  radius  values  approaching  Tmaxi  and  rmax2  is  unchanged  from  the  terrestrial 
interference  case.  The  performance  of  the  ML  approach  was  significantly  better  than  the 
other  two  approaches  and  remained  relatively  constant  for  all  array  sizes. 
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(a) 


(b) 

Figure  26:  Effect  of  external  noise  on  the  Watson-Watt,  Butler  matrix,  and  ML  ap¬ 
proaches  using  (a)  a  four-element  and  (b)  an  eight-element  circular  array  with  a  radius 
r  =  0.2 A  and  a  signal  bearing  of  0  —  =  0°. 
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(a) 


(b) 

Figure  27:  Effect  of  the  antenna  array  radius  on  the  performance  of  the  Watson- Watt, 
Butler  matrix,  and  ML  approaches  using  an  eight-element  circular  antenna  array  in  the 
presence  of  external  noise  {snr  =  10  dB)  using  (a)  a  four-element  circular  array  with  a 
signal  bearing  of  <f)  — 'd  =  45°  and  (b)  an  eight-element  circular  array  with  a  signal  bearing 
of  ^  -  t?  =  0°. 
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Figure  28:  Effect  of  external  isotropic  noise  on  the  Watson- Watt,  Butler  matrix,  and 
ML  approaches  using  an  eight-element  circular  antenna  array  with  a  radius  r  =  0.2A. 


Figure  29:  Effect  of  the  antenna  array  radius  on  the  performance  of  the  Watson- Watt, 
Butler  matrix,  and  ML  approaches  using  an  eight-element  circular  antenna  array  in  the 
presence  of  external  isotropic  noise  {snr  =10  dB). 
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4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  report,  the  Watson- Watt  approach  and  Butler  matrix  approach  -  two  direction 
finding  methods  developed  for  four  and  eight-element  circular  antenna  arrays  and  pro¬ 
posed  for  wideband  applications  -  were  evaluated  theoretically  and  through  simulations 
for  error  causing  conditions  found  in  the  real  world.  These  conditions  include  internal  re¬ 
ceiver  noise,  external  terrestrial  and  isotropic  noise,  cochannel  interference  and  multipath. 
Additionally  comparisons  through  simulation  were  also  carried  out  using  a  maximum  like¬ 
lihood  (ML)  approach  optimized  for  a  single  signal  in  order  to  determine  the  performance 
trade-off  incurred  using  either  of  the  two  proposed  wideband  DF  techniques.  The  ML  ap¬ 
proach  is  most  closely  related  to  correlative  interferometer  or  vector  matching  approaches 
used  in  narrowband  DF  systems  (e.g.  the  Watkins-Johnson  8996-1  DF  processor  or  the 
Rohde  &  Schwarz  DDF  190),  the  main  difference  being  that  these  approaches  do  not 
compensate  for  the  statistical  distribution  of  noise  whereas  ML  does. 

From  the  implementation  standpoint,  the  Butler  matrix  approach  is  preferable  since 
it  requires  only  two  receiver  channels  compared  to  three  receiver  channels  for  the  Watson- 
Watt  approach.  This  reduces  the  expense  in  terms  of  wideband  receivers,  digitizers, 
etc.  although  some  additional  circuitry  in  the  front-end  is  required  including  a  quadra¬ 
ture  mixer.  By  way  of  comparison,  the  ML  approach  would  require  a  receiver  for  every 
antenna  in  the  array  plus  a  more  sophisticated  bearing  processor  owing  to  the  greater 
computational  complexity  of  this  approach. 

In  terms  of  the  effects  of  internal  receiver  noise,  both  the  Watson- Watt  and  Butler 
matrix  approaches  had  similar  performance  and  actually  outperformed  the  ML  approach 
for  smaller  antenna  arrays  (i.e.  less  than  one  fifth  of  a  wavelength)  in  simulations.  This 
surprising  result  was  due  to  the  fact  that  the  ML  approach  uses  more  receivers  and  is 
therefore  exposed  to  a  greater  amount  of  noise  as  a  result.  The  simulations  also  highlighted 
problems  with  certain  antenna  array  sizes  where  channel  gains  of  some  of  the  channels 
go  to  zero  resulting  in  gross  bearing  errors.  This  limits  the  maximum  antenna  array  size 
that  can  be  used,  or  equivalently,  the  upper  frequency  of  the  array.  There  is  no  such  gain 
problem  for  the  ML  approach,  although  the  four-element  array  suffers  ambiguity  problems 
for  the  same  radii;  the  eight-element  array  was  unaffected  however. 

In  terms  of  the  effect  of  external  influences,  including  terrestrial  or  isotropic  noise, 
cochannel  interference,  and  multipath,  the  Watson- Watt  approach  was  twice  as  accurate 
as  the  Butler  matrix  method.  The  main  reason  for  the  poorer  accuracy  of  the  Butler 


56 


matrix  method  is  that  errors  in  the  reference  channel,  which  is  required  for  quadra¬ 
ture  demodulation  of  the  directional  channel  into  its  sin  <f>  and  cos  (j)  bearing  components, 
translate  directly  into  bearing  errors.  In  the  Watson- Watt  method,  the  sin  (f>  and  cos  <f> 
bearing  components  are  contained  in  separate  channels  so  that  the  reference  channel  is 
only  required  to  resolve  180°  bearing  ambiguities. 

Generally,  the  effect  of  external  influences  on  accuracy  was  constant  with  increas¬ 
ing  antenna  array  size  or  accuracy  got  worse.  For  the  Watson- Watt  and  Butler  matrix 
approaches  there  also  appeared  to  be  very  little  performance  advantage  in  using  an  eight- 
element  antenna  array  compared  to  a  four-element  antenna  array.  The  problem  of  critical 
array  sizes  resulting  in  zero  channel  gain  and  causing  gross  bearing  errors  was  also  found 
to  limit  the  usable  size  of  the  antenna  array/frequency  range.  The  ML  approach  was 
found  to  provide  significantly  better  performance,  particularly  when  using  the  eight  ele¬ 
ment  array,  indicating  that  performance  is  related  to  the  number  of  receiver  channels  used 
(i.e.  the  Butler  matrix  approach  which  used  the  fewest  also  performed  the  most  poorly, 
whereas  the  ML  approach  which  used  the  most  performed  the  best). 

For  most  practical  applications,  external  influences  are  more  likely  to  play  a  greater 
role  in  degrading  DF  performance  than  internal  noise.  Hence  performance  under  these 
conditions  is  more  critical  in  assessing  the  overall  performance  than  internal  noise  effects 
(although  internal  noise  effects  cannot  be  ignored). 

The  issue  of  critical  antenna  array  sizes  can  be  broken  into  two  problems:  zero  gain 
in  the  reference  channel  for  an  circular  array  radius  of  approximately  0.3  —  0.4A,  and  zero 
gain  in  the  directional  channels  for  a  circular  array  radius  of  approximately  0.5  —  0.6A 
(where  A  is  the  signal  wavelength  at  the  desired  maximum  frequency).  The  solution  to  the 
reference  channel  zero  gain  is  to  use  an  extra  antenna  at  the  center  of  the  antenna  array 
rather  than  from  the  summation  of  the  circular  array  antenna  outputs.  This  could  be 
problematic,  however,  if  the  central  position  is  occupied  by  the  supporting  mast.  There 
is  no  obvious  way  to  deal  with  the  zero  gain  problem  in  the  directional  channels  so  that 
the  upper  limit  on  circular  array  size  is  approximately  0.5  —  0.6A.  The  ML  approach,  by 
comparison,  performs  well  for  circular  array  radii  up  to  0.7A  without  requiring  an  extra 
central  antenna. 

Based  on  these  findings,  the  main  trade-off  is  accuracy  versus  cost  and  speed.  The 
ML  approach  achieves  the  highest  accuracy  under  most  conditions  and  has  the  greatest 
frequency  range  of  operation.  It  also  requires  the  most  number  of  receiver  channels  (as 
many  as  the  number  of  antennas)  and  has  the  greatest  processing  requirements  making 
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realtime  operation  a  much  more  expensive  proposition.  The  Watson- Watt  approach  is 
less  accurate,  but  requires  only  three  receiver  channels  and  has  very  simple  processing 
requirements  making  it  a  much  more  cost  effective  approach.  The  Butler  matrix  approach 
reduces  the  number  of  receiver  channels  even  further  to  two  and  also  has  simple  processing 
requirements,  but  sacrifices  more  accuracy,  having  RMS  bearing  errors  twice  those  of  the 
Watson- Watt  approach  under  practical  conditions. 

The  conclusion  is  that  for  a  system  where  not  only  accuracy,  but  also  low  cost,  and 
realtime  operation  are  the  major  requirements,  the  Watson- Watt  approach  is  the  best 
choice  for  current  applications.  In  the  future,  however,  this  choice  may  need  to  be  re¬ 
assessed  as  lower  cost  receivers  and  faster  processing  equipment  make  the  ML  approach 
a  more  attractive  alternative. 
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